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' Abstract 
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We use spin-density functional theory to obtain novel expressions for the charge and spin local-field factors of an 
^ 'I electron gas in terms of its electron-pair structure factors. These expressions (i) satisfy the compressibility and spin 
susceptibility sum rules; (ii) keep account of kinetic correlations by means of an integration over the coupling strength; 
and (iii) provide a practical self-consistent scheme for evaluating linear response and liquid structure. Numerical 
illustrations are given for the dielectric response of the paramagnetic electron gas in both three and two dimensions. 
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PACS numbers: 71.10.Ca - Electron gas 

Key Words: C. Electron-electron interactions; C. Dielectric response; C. Spin dynamics 



An interacting electron gas on a uniform neutralizing background is the reference system for many calculations 
of electronic structure in condensed matter systems [1]. Its dielectric and spin susceptibilities provide key inputs in 
density functional theory [2] and in studies of quasi-particle properties such as the effective mass and the effective 
Lande g- factor [3]. These properties are now known from experiment for carriers in semiconductor structures over a 
wide range of density [4]. A great deal of relevant information has been coming from Quantum Monte Carlo (QMC) 
studies [5-9], but theoretical understanding continues to attract interest. 
ON ' The issue of thermodynamic consistence in the theory of linear response for an electron gas, though raised in the 
early work of Vashishta and Singwi [10] through an ad hoc scheme to account for the compressibility sum rule, has 
ON ! not received much attention. In brief, thermodynamic consistence requires that the linear susceptibilities (at zero 
temperature) should reduce in the static, long-wavelength limit to the values obtained from the second derivatives of 
| the (ground-state) energy. Of course, it is possible to model the susceptibilities by parametrized analytical expressions 
that embody these limiting behaviors as they are known from QMC data [11, 12]. In this Letter we propose instead 
a fully self-consistent theoretical scheme, which uses the mean particle and spin densities of the system as the only 
inputs for the calculation of the susceptibilities and of the electron-pair structure. Our approach is aimed at satisfying 
both the compressibility and the spin susceptibility sum rule: indeed, the QMC studies show that spin correlations 
and spin ordering are enhanced on increasing the coupling strength. Our starting point is the expression of the 
exchange-correlation energy E KC [n^, n\\ as a functional of the spin densities n^(r) and «j.(r) in an inhomogeneous 
d-dimensional electron gas with interaction potential v(r) = e 2 /r, 



E ^ n h n l] = J^2 J q dX Y.J ddr J d d r'v(\r-r'\) [ {V x \8n a (r)5n a , (r')\V x ) - 8 aa ,n a (v) 5(r - r') 



(1) 



where \^B\) is the ground state of the gas with interaction potential v\(r) = X/r and 5ri a (r) = (r) — n a (r) . The spin- 
density operators n a (r) are given by n a (v) = ipl-(r)tfj a (r) in terms of Schrodinger field operators obeying canonical 
anticommutation relations. In deriving Eq. (1), which is obtained by an immediate extension of the treatment given by 
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Dreizler and Gross [2] for the paramagnetic state, the technique of coupling-constant integration has been adapted to 
the context of density functional theory in order to account for kinetic correlations. Upon introducing as independent 
variables the total particle density n(r) = n-f(r) + n^(r) and the spin magnetization density m(r) = n-f(r) — nj,(r), 
Eq. (1) can be written in the more convenient form 



E* c [n,m] = \ J d d v J d d r' v(\r - r'|) n(r)n(r') [g(r, r') - 1] 
where the coupling-averaged pair function g(r,r') is defined by 

i r 2 

s(r,r') = — / dX g\(r,r'). 
e Jo 



(2) 



(3) 



Here, g\(r, r') is the pair distribution function in the gas with interaction potential v\(r) of strength A and is, of 
course, a functional of n and to. 

We are in fact interested in the linear-response regime, where the inhomogeneity of the system is due to a weak 
static perturbation. We can then introduce the exchange-correlation kernels for the homogeneous electron gas as 



K™\n,m;\r-v'\ 



6 2 E X 



(5n (T (r)(5n (T '(r') 



hg 



and the local field factors G aC r> (q) as 



Vq G a 



'(<?) 



d d vK™ (r) exp(-iq-r) , 



(4) 



(5) 



with v q the Fourier transform of the electron-electron potential. In Eq. (4), after taking the second functional 
derivative the inhomogeneous spin densities are replaced by those of the homogeneous gas, corresponding to a particle 
density n and a magnetization density to. 

While the formulation given so far is quite general (within neglect of spin-orbit coupling) , in the following we restrict 
ourselves for simplicity to the paramagnetic state (i.e. fh = 0). In this case the matrices in Eqs. (4) and (5) can be 
diagonalized and expressed in terms of the charge and longitudinal-spin response kernels, 



and 



K+(n;\r-r'\) = 



^xct™; |r- r' 



S 2 E X 



Sn(r)Sn(r') 



S 2 E X 



5m(v)5m(v') 



hpg 



hpg 



(6) 



(7) 



The corresponding local field factors are G+(q) and G-(q), say. Notice that in the calculation of K~ c (n;r) the 
functional derivatives in Eq. (7) must be evaluated on the spin-polarized gas before taking the paramagnetic limit. 

The evaluation of the functional derivatives in Eqs. (6) and (7) can be explicitly carried out from Eqs. (1) and 
(2). The compressibility and spin susceptibility sum rules can be satisfied by making at this point a long- wavelength 
approximation. This involves, as an instance, the replacement 5g(ri, r')/Sn(r) — > S(ri — r)dg(\r — r'\)/dn. Some 
straightforward algebra leads to the expressions 



g(r) - 1 + 2h 



dg(r) 
dn 



1 - 2 d 2 9(r) 



+ 2 U 



dn 2 



and 



K xc (n;r) = — 



2 n ' 



d 2 g(r) 



dm 2 



m=0 



(8) 



(9) 
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The local field factors are thus given in this approximation by 

d_ 1 

dn 2 dn 2 



G + (?)=(l + 2n^ + in 2 ^)G( 9 ) ,10, 



and 



where we have defined 



(ii) 



m=0 



G(q) 



■^l' dX I^ v i ls ^-^- 1] - (12 » 
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Here, Sa(<z) is the static structure factor of the electron gas with a pair potential v\{q) given by AirX/q 2 in d = 3 and 
by 2ttA/(7 in d = 2. That is, 

S x (q) = l + nf rf d r[ ffA (r)-l]cxp(-zq.r) . (13) 

Eqs. (10) - (12) arc the main new results of the present work. 

Let us pause at this point to explicitly show that Eqs. (10) - (13) satisfy the two thermodynamic sum rules. By 
taking the long-wavelength limit q — > inside the integral in Eq. (12) we find 

Um G() = .MM 

where 

e xc (n,m) = ± £ dX J 0- d v q , [S x (q') 1] (15) 

is the exchange-correlation energy in the homogeneous, spin-polarized electron gas. Since e xc is equal to the difference 
e gs — £o between the ground-state energies of the real and of the non-interacting electron gas, it follows from Eqs. (10) 
and (14) that 

!i™ \G + (q) = J- ( -L - I) (16) 



q-*o n 2 v q \K K 

K and K being the corresponding compressibilities. Similarly, from Eqs. (11) and (14) we get 



lim G-(q) = — 



1 d 2 s xc (n,0 



q^O nv q d£ 2 

with £ = fh/n and hence 



(17) 
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limG_(g) = ^ (18) 

9^0 Vg \XP XsJ 

where [Ib is the Bohr magneton, \p the Pauli susceptibility and Xs the spin susceptibility of the interacting electron 
gas. Equations (16) and (18) are the correct relations between the long-wavelength limit of the static response of the 
paramagnetic electron gas and the second derivatives of its ground state energy [12]. 

Returning to Eqs. (10) - (12), they evidently need supplementing by a method for the evaluation of the electron- 
pair structure. Such a calculation can be carried out self-consistently by one of two alternative routes: (i) by using 
the static linear response functions to construct effective spin-dependent electron-electron interactions that yield the 
electron-pair distribution functions through the solution of an electron-electron scattering problem [13 - 15]; or (ii) by 
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constructing dynamic susceptibilities from the static local-field factors and imposing that the fluctuation-dissipation 
formula relating the electron-pair distribution function to the van Hove dynamic structure factor is satisfied [16]. The 
former method has already been tested to some extent and shown to yield accurate results for quantum plasmas over 
a physically significant range of coupling strength [17]. In the following we present and discuss the second method 
mentioned above, in regard to the dielectric response and charge-charge correlations for electrons in d — 3 and d = 2. 
The study of spin-spin correlations requires entering the spin-polarized state, as already mentioned, and is left for 
future work. 

In brief, following the work of Singwi et al. [16] we write an approximate expression for the dynamic dielectric 
susceptibility Xc(l, from the static local-field factor G+(q) as 

Xc(g,w) = - r — , — -. r (19) 

l-v q [l-G + (q)]xo{q,u)) 

where Xo(l,u) is the dynamic (Lindhard-Stcrn) susceptibility of the ideal Fermi gas. The electron-pair structure 
factor can then be calculated as 

S{q) = ~ f dculmxdq^). (20) 

Self-consistence between pair structure and dielectric response is ensured by the fact that the local-field factor to be 
used in Eq. (19) depends on the pair structure according to Eqs. (10) and (12). Of course, Eqs. (19) and (20) hold 
for each value of A and must be solved together with Eqs. (10) and (12) up to the physical value of the coupling 
strength in the system. 

The results of this calculation for a weakly coupled electron gas in d = 3 arc shown in Figures 1 and 2. It is seen 
from Figure 1 that the local-field factor G+(q) obtained from the present theoretical approach is quite close to the 
QMC data of Moroni et al. [7] over a substantial range of values for q/JtF- Important deviations from these data 
are found only at high momenta q > 2kp, where the details of the dielectric response have little relevance for most 
practical purposes: in fact, in our approach the asymptotic value of G(q) at large momenta is a constant determined 
by the quantity 1 — g(r = 0), rather than showing a leading term of order q 2 . It is also seen from Figure 1 that 
our self-consistent results are in good agreement with the parametrized expression of G + (q) given by Ichimaru and 
Utsumi (IU) [18], who combined QMC data on the compressibility with a ladder diagram calculation of g(r = 0) and 
a previous dielectric formulation by Utsumi and Ichimaru [11]. The peak of G+(q) for momenta approaching 2kp, 
which in their formulation arose from a logarithmic term associated with exchange, comes in our approach from a 
logarithmic factor remaining inside the integral in Eq. (12) after angular integration in d = 3. This peak is absent in 
the so-called STLS dielectric theory of Singwi et al. [16], as can again be seen from Figure 1. 

Our results for the pair distribution function g(r) at full self-consistence in d = 3 are reported in Figure 2, in 
comparison with QMC data kindly provided by P. Ballone (unpublished) and with the results of the STLS theory. 
Evidently, the present results for the electron-pair structure are excellent in this weak-coupling regime. With increasing 
coupling strength the above-mentioned peak in G+(g), however, introduces unwanted oscillations in g{r). It appears 
that as the coupling increases the frequency dependence of the local-field factors [19] may become important in the 
calculation of pair correlations via the fluctuation-dissipation theorem. Much more satisfactory results are instead 
obtained by calculating these correlations through the solution of an electron-electron scattering problem with effective 
scattering potentials embodying the many-body effects from the static local-field factors [17]. 

Finally, Figure 3 reports similar results for the local-field factor G+(q) in a weakly coupled electron gas in d = 2. In 
this case satisfying the compressibility sum rule ensures that the local-field factor from QMC data [12] is reproduced 
almost up to q ~ 2kp. Again a peak is found to be present in G + (q) for q ~ 2kp, which is absent in an STLS 
formulation. 

In summary, we have proposed a self-consistent theoretical approach to the evaluation of the local-field factors 
entering the dielectric and longitudinal-spin response of an interacting electron gas and given some examples of 
the results that it predicts for the dielectric response and the charge-charge correlations at weak coupling. In the 
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present formulation the form of the local-field factors is very accurate in the most important range of momenta where 
the thermodynamic sum rules dominate the physical behavior, but may still be too sensitive to exchange effects at 
momenta of order 2kp. For this reason, structural self-consistence appears to be best enforced at strong coupling 
through the use of the local-field factors in constructing effective many-body scattering potentials. 
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FIG. 2. The radial distribution function g(r) of a paramagnetic electron gas in d = 3 at coupling strength r a = 1, as a 
function of reduced distance rkF. 
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